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We consider the ordering kinetics of a nonconserved scalar field advected by a uniform shear flow. 
Using the Ohta-Jasnow-Kawasaki approximation, modified to allow for shear-induced anisotropy, 
we calculate the asymptotic time dependence of the characteristic length scales, Ln and L±, that 
describe the growth of order parallel and perpendicular to the mean domain orientation. In space 
dimension d = 3 we find Lu ~ 7i 3//2 , L± ~ t 1//2 , where 7 is the shear rate, while for d = 2 we find 
JLrii ~ 7 1/ ' 2 t(ln£) 1 ' 4 , L± ~ 7~ 1,/2 (lnf) _1 / 4 . Our predictions for d — 2 can be tested by experiments 
on twisted nematic liquid crystals. 



The coarsening dynamics of systems quenched from 
a disordered phase into a two-phase region are by now 
reasonably well understood JO. Domains of the ordered 
phases form rapidly, and then coarsen with time t, i.e. 
there is a characteristic length scale ('domain size'), L(t), 
which grows with time, typically as a power law. Fur- 
thermore, there is good evidence for a form of dynamical 
scaling in which the domain morphology is statistically 
scale- invariant when all lengths are scaled by L(t) JjJ. 

In recent years attention has been directed at sys- 
tems subjected to external driving, e.g. an imposed shear 
. The shear induces an anisotropy, leading to differ- 
ent coarsening rates in directions parallel and perpendic- 
ular to the flow. This has been observed in experiments 

An important open question for such driven systems is 
whether the coarsening continues indefinitely (for an infi- 
nite system), as in the case of no driving, or whether the 
driving force arrests the coarsening, leading to a steady 
state. For the case of a sheared phase-separating binary 
fluid, it has been argued, on the basis of the stability of a 
single drop of one fluid immersed in another, that the do- 
main scale will eventually saturate at a maximum length 
scale, L max , determined by the shear rate: L max ~ cr/777, 
where a and 77 are the surface tension and viscosity re- 
spectively 0. In the multi-domain context, it has been 
argued H that a steady state would be achieved through 
the shear-induced stretching and breaking of domains. 
However, the experimental evidence for a steady state is 
not completely clear-cut. In particular, as emphasized by 
Cates et al. H, saturation of the domain length occurs 
naturally in a finite-size system when the domain length 
becomes of the same order as the system size. 

A second important question concerns the nature of 
the growth laws, and the nature of scaling - if it exists 
- in this anisotropic system. Naively, one might expect 
two characteristic length scales, L11 and L±, measuring 
correlations parallel and perpendicular to the flow. In 
dimension d = 2, for example, it would be natural to 
conjecture that the coarsening follows conventional scal- 
ing when lengths along and perpendicular to the flow 
are scaled by these characteristic lengths. We will show, 
however, that the actual scaling is subtly different. 



To address these issues it would be helpful to have an 
exactly soluble model. The large-n limit of the n-vector 
model has been solved for a conserved order-parameter 
field, but without hydrodynamics (i.e. the order param- 
eter is simply advected by the shear flow) Rl. Length 
scales L11 and Lj_, growing as j(t 5 / lni) 1 / 4 and (t/ lnf) 1 / 4 
respectively, describe correlations along and perpendicu- 
lar to the flow. Scaling is not strictly satisfied (there is 
instead a form of 'multiscaling'), but this is presumably 
an artifact of the large-n limit, as in the zero-shear case 
[^1 . There is no saturation at late times in this model. 
However, for the large-n model there are no domains, so 
the concept of stretching and breaking loses its meaning. 
What is needed is an analytically tractable model for a 
scalar order parameter. For a nonconserved field, the 
Ohta-Jasnow-Kawasaki (OJK) approximation fulfills 
this requirement: In the absence of shear, it captures the 
essential features of the coarsening process [jlj . 

In this Letter we present the results of applying the 
OJK approach to a nonconserved scalar field advected 
by a uniform shear flow. The equation of motion for the 
order parameter is 

+ V •(!*£) = V 2 cj)~V'(cb), (1) 

where V(<f>) is a symmetric double well potential, and 
u = 7yx is the velocity of the imposed shear flow, with 
the flow in the x-direction. Two aspects of real binary 
fluids are neglected in this model: the order parameter 
is not conserved, and it is simply advected by the shear 
rather than being coupled to the fluid velocity through 
the Navier-Stokes equation. This approach is, however, a 
very instructive first step which yields important insights 
concerning the main questions raised earlier, namely the 
nature the asymptotics and of dynamical scaling, in the 
physically correct context of a scalar field. Furthermore, 
in d = 2 this model describes the coarsening dynam- 
ics of a twisted nematic liquid crystal, in which discli- 
nation lines, separating domains of opposite twist, relax 
viscously, driven by their line tension p0| . Under shear, 
this system will furnish an experimental test of our pre- 
dictions. It should be noted that our analysis leads to 
very different behavior in d = 2 and d — 3: This is one 
of the main results of the present work. 
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The result of the OJK analysis is that, for space di- 
mension d — 3, Lii ~ 7 f 3 / 2 and L±_ ~ t 1 / 2 , i.e. the coars- 
ening rate parallel to the flow is enhanced by a factor 
"ft relative to the unsheared system, while the growth of 
L± is unchanged (the same features describe the large-n 
result Q). Hence there is no saturation of the coarsen- 
ing. For d = 3, furthermore, conventional scaling holds. 
For d = 2, however, very different results are obtained: 
L|| ~ 7 1 / 2 £(lni) 1 / 4 , and L± ~ 7 - 1 / 2 (ln<)- 1 / 4 . These re- 
sults imply L\\L± ~ i for d = 2, i.e. the product of the 
two length scales, or 'scale area', is independent of the 
shear rate and has the same form as for the unsheared 
system, where Lii = Lj_ ~ t 1 / 2 . We will show that this 
result can be understood by a topological argument. An 
important subtlety in d = 2 is that Ln and L± have to be 
defined as characteristic scales parallel and perpendicu- 
lar to the mean domain orientation, instead of the flow 
direction. This distinction is not important for d = 3, 
but crucial for d = 2. Only with the new definition is 
dynamical scaling recovered. For d = 2, furthermore, L±_ 
decreases with time asymptotically. Our approach breaks 
down, however, when L± becomes comparable with the 
width, £, of the interfaces between domains, which occurs 
after a time of order exp(const/7 2 £ 4 ), when we expect 
the domains to break, as observed in simulations with 
conserved dynamics in d = 2 [El . 

The OJK approach starts from the Allen-Cahn equa- 
tion pH relating the normal component of the inter- 
face velocity, v n , to the local curvature of the interface, 
K = V • n, where n is the normal to the interface, 



d t m + jyd x m = V 2 m - D ab (t)d a d b m 



-V • n + u • n, 



(2) 



where the final term is the drift due to the shear. The 
derivation of this equation from (|]) follows the same 
route as the zero-shear case The next step is to 
introduce a smooth auxiliary field m(x, t) whose zeroes 
coincide with those of <j>. In a frame comoving with the 
interface one has dm/dt = = dtm + v n \Vm\. Combin- 
ing this with (0), and using 



n = Vm/|Vm|, 

yields the following equation for m, 

d 

dtm + u ■ Vm = V 2 m — n a n b d a d b m 

a, 6=1 



(3) 



(4) 



So far this is exact. Equation (0) is highly non-linear, 
however, due to the implicit dependence of n on m 
through (^). The OJK approximation involves lineariz- 
ing the m equation by replacing the product n a rib by its 
spatial average, 

D a b{t) = (n a (x,t)nb(x,t)) , (5) 

leading to the following equations for m and D: 



D ab = 



d a m d b m 
(Vm) 2 



(6) 
(7) 



In the absence of shear the coarsening is isotropic, and 
the matrix D has the simple form D ab = S ab /d, indepen- 
dent of t. The resulting diffusion equation for m is readily 
solved, leading to a \ft coarsening. For 7 7^ 0, the coars- 
ening is anisotropic: D ab is both non-diagonal and time- 
dependent. From the definition (||) of D ab , though, the 
sum- rule, Tr_D = 1, is trivially valid. From the symme- 
try of (^j|) under the combined transformations x — » —x, 
y — > — y at fixed z, and under the separate transformation 
z — > — z at fixed x, y, we see that D ab has a block diago- 
nal form, with D xy = D yx the only non-zero off-diagonal 
elements. 

In Fourier space, (^) reads 

^^-7k x ^^ = -^n ab (t)k a k b rn(Kt), (8) 



ab 



Clab{t) = 5 ab - D ab (t) 



(9) 



This is readily solved by the change of variables q = Ak, 
t = t, and /i(q, r) = m(k, t), where A has elements 



A a b = 5ab + lt8 a2 8 b i. 



(10) 



In the new variables, the left-hand side of (g) becomes 
<9 r /i(q, t). Integrating the equation, and transforming 
back to the original variables, gives 

m(k,t) =m(k(*),0)exp ^-i ^ k a M ab (t) k b ^j , (11) 

where k(t) = (k x ,k y + r yk x t,k z ) and the matrix M is 
given by 



M(t) = A T (t) R(t) A(t) 
R{t)=4 I dt'[A T {t')]- 1 rt(t')[A(t')]- 1 . 



(12) 
(13) 



Equations (|ll|-|l3|) determine the function m(k, t) com- 
pletely if the matrix D ab {t) is known. However, D is itself 
determined from the distribution for m, via (]?]), so we 
have to solve these equations self-consistently. We take 
the initial condition, m(k, 0) to be a gaussian random 
variable with correlator (m(k, 0)m(k', 0)) = A <5(k + k'). 
Then the real-space correlation function of m, G(r,t) = 
(m(x + r, t)m(x, t)) is obtained from (11) as 



G(r,i) = G(0,t) exp ^(M" 1 )^ r b j 



(14) 



where the precise expression for G(0, t) is not relevant 
in what follows. All the information concerning domain 
growth in this system is contained in the matrix M ab (t). 
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To obtain a closed set of equations we have to express 
the elements of the matrix D in terms of the elements 
of M. To to this we first rewrite ([?]) using an integral 
representation for the denominator: 



D„ 



du l^d a mdbm exp ^— ^(Vm) 2 Jj . (15) 



Since m is a Gaussian field, the required average can 
be computed using the probability distribution -P(v) of 
the vector v = Vm. This distribution is determined 
by the correlator {v a v b ) — — d a df, (m(r, t) m{0,t))\ T =o = 
(rn 2 ) (M ~ 1 ) a b, from which one infers that 



P(v) oc exp 



1 



2(m 2 



v * M <> 



bV b 



(16) 



Carrying out the average in ( |l5| ) gives 



Dab 
N ab 



, /detAf\ 1/2 , AT u 



1 

/ du . 

2 J \detN 

M ab + U Sab 



(17) 
(18) 



We now proceed to outline the solution of this closed 
set of equations. Full details will be given elsewhere Jl2| . 
We are interested in the large-i asymptotics. This limit 
simplifies the analysis which is still, however, quite sub- 
tle. The results are very different in three and two di- 
mensions, so we discuss these cases separately. 

We begin with some general remarks on the expected 
form of D a b- The effect of the shear is to produce elon- 
gated domain structures aligned, at late times, at a small 
angle (~ 1 /jt - see below) to the flow direction. As a re- 
sult, the component n\ of the normal to the interface 
is very small almost everywhere at late times, implying 
Dn — > for t — ► oo. The sum rule, TrD = 1, is there- 
fore exhausted by the remaining diagonal components of 
D for large t. In particular, for d = 2 we have D 2 2 — > 1 
for t — ► oo, and fin — * in (g). This case requires special 
care. For d = 3, on the other hand, D22 + -D33 — * 1, and 
it turns out that both D22 and D33 approach non-zero 
limits. This case is, therefore, simpler to analyse. 

The case d = 3. With the assumption that Du — ► 
and D12 ~ 1/t, while D22 and D33 remain non-zero for 
t — ► 00, the asymptotics of the matrix elements M ab are 
readily obtained from (§),©, @ and @: 



7^(1 



M 12 = 2 7 t 2 (1 - Df 2 ) 
M 22 = 4t(l-Dg) 
M 33 = UD% , 



(19) 



where D^ is the large- i limit of D22, while M13 = M23 
by symmetry. 



Using these limiting forms, the integrals ([i~7|) can be 
evaluated asymptotically. After some algebra one finds 



(20) 




while = 1- 



and D13 = D23 = 0. Equation 
(|20|) has the non-trivial solution DS^ — 4/5, implying 
^33 = 1/5 and demonstrating the self-consistency of the 
initial ansatz. Finally we find Du — > 31n(7i)/(7i) 2 and 
Z?i2 — > —6/(57*), consistent with our initial assumption. 

The characteristic length scales in the system are given, 
from ( |l4| ) , by the square roots of the eigenvalues of the 
matrix M. Using |[|), with D%% = 4/5, we find 



L« = 



It 



3/2 



T - 1 f 1 / 2 



1/2 



for t — > 00. The corresponding eigenvectors are 




(21) 



(22) 



implying that the principal axes in the xy plane are ro- 
tated, relative to the x and y axes, through an angle 
3/27*, which we can interpret as the angle between the 
average orientation of the domain structure and the flow 
direction. 

The three length scales are all distinct, though L± and 
L3 grow in the same way, and coarsening continues indef- 
initely - the system does not approach a stationary state. 
The matrix elements M a b grow in the way expected if a 
naive form of scaling holds: Mu = L 2 , M\2 ~ L\\L±, 



M22 ~ and M33 = L3. This means that dynamical 



scaling holds, and the scaling variables can be taken to 
be x/Lu, y/Lj_, and z/L^. The same simple structure 
does not, however, hold in d = 2. 

The case d = 2. For d = 2 the self-consistency prob- 
lem is more tricky, because the quantity 1 — £*22^) = 
-Dn(i) tends to zero as t — > 00. For d = 2 the inte- 
grals ([It]) can be evaluated exactly. Making the assump- 
tions, to be verified subsequently, that asymptotically 
Mu > Mvi > M 22 , that TrM > Vdet M > Af 22 , 
and that D\2 ~ 1/t, one can derive the following self- 
consistent equation for Du{t)'- 



1 ( gdt^Du^) 
Dn{t) -^[l^Du(t') 



1/2 



with asymptotic solution 
Bii(i) 



27t v / In7l' 



(23) 



(24) 



Using this result, the asymptotic results for the matrix 
elements, and the determinant, of the correlation matrix 
M are obtained as 
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Mn(t) = A-yt 2 ^ln-ft 
M 12 (t) = At^/kryt 
M 22 (t) = (4/ 7 ) y/k^~t 
detM(i) = 4i 2 , (25) 

while D12 = —I /"ft. These results confirm, a posteriori, 
the assumptions made in their derivation, i.e. the solu- 
tion is self-consistent. The asymptotic results for the M a & 
seem to imply detM = 0, in contradiction to (p5[). To 
obtain (^5|) one has to keep subdominant contributions 
to the M ab . [Note that detM = deti? from @]. 

The characteristic length scales are given, as before, 
by the square roots of the eigenvalues of M. The eigen- 
values of any 2x2 matrix can be expressed as X± = 
[Tr M ± ,/(Tr M) 2 — 4 det M}/2. Using Tr M > Vdet M 
we obtain A + = Tr M and A_ = detM/(TrM), whence 

L,| = 2^t (ln 7 ty/*, L ± = 1 . (26) 

1 /7(ln7<) l /4 

The corresponding eigenvectors are 

*« = (!)' —{-?)• < 27 » 

giving a tilt angle 1/7^ between the domain orientation 
and the flow direction. 

The scale area in two dimensions is LuL± = 2t, inde- 
pendent of 7. This is the same result as the zero-shear 
case, where D a b — 5 a b/d implies M a b = 2tS a b for d = 2, 
i.e. L(t) = \/2~i. This result is special to d — 2 and can be 
understood as follows. For an isolated domain, the rate 
of change of the area enclosed by the domain boundary is 
dA/dt = §dlv n = § dl (u - V) n from (§). The second 
term is a topological invariant, equal to — 2n from the 
Gauss-Bonnet Theorem, while the first term is equal to 
f. d 2 x V • u, which vanishes for any divergence-free shear 
flow. While the 7-independence of dA/dt has been proved 
only for closed loops of domain wall, we expect a similar 
result to hold for the scale area, i.e. d(L\\L±)/dt = const. 
It is very encouraging that the OJK approximation cap- 
tures this feature of the d — 2 coarsening. There is no 
equivalent result in d = 3 because the surface integral 
/ dS V • n is no longer a topological invariant. 

The scaling is nontrivial in d = 2 because, although 
Mn — L 2 , consistent with naive scaling, the correspond- 
ing results M12 ~ L\\L±, M 2 2 ~ L\, found in d = 3 no 
longer hold in d = 2. This is associated with the fact 
that the leading-order contribution to detM vanishes. 
The consequence is that scaling only holds when referred 
to the unique scaling axes (|27]), which are themselves 
time-dependent. 

The most interesting and suggestive feature of the 
d = 2 result is that L± tends to zero as t — > 00. Since our 



treatment is based on the 'thin wall' limit, in which do- 
main walls are treated as having zero width, it will break 
down when L± becomes comparable with the width, £, 
of the walls, at which point we conjecture that domains 
will break, possibly arresting the coarsening. This can be 
tested by experiments on twisted nematic liquid crystals. 
In d — 3, the present work provides strong evidence that, 
at least for the nonconserved scalar field considered here, 
the coarsening state proceeds indefinitely. In this respect 
it is interesting that, in their d = 3 simulations (including 
both conservation of the order parameter and hydrody- 
namics), Cates et al. H found no evidence for a steady 
state structure emerging that is independent of the sys- 
tem size, i.e. observed steady states could be attributed 
to finite size effects. 

We thank Peter Sollich for a useful discussion. This 
work was supported by EPSRC (UK) under grant 
GR/L97698. 



[1] A. J. Bray, Adv. Phys. 43, 357 (1994), and references 
therein. 

[2] A. Omiki, J. Phys.: Condens. Matter 9, 6119 (1997), and 
references therein. 

[3] C. K. Chan, F. Perrot, and D. Beysens, Phys. Rev. A 43, 
1826 (1991); A. H. Krall, J. V. Sengers, and K. Hamano, 
Phys. Rev. Lett. 69, 1963 (1992); T. Hashimoto, K. Mat- 
suzaka, E. Moses, and A. Onuki, Phys. Rev. Lett. 74, 126 
(1995); J. Lauger, C. Laubner, and W. Gronski, Phys. 
Rev. Lett. 75, 3576 (1995). 

[4] D. H. Rothman, Phys. Rev. Lett. 65, 3305 (1990); P. 
Padilla and S. Toxvaerd, J. Chem. Phys. 106, 2342 
(1997); A. J. Wagner and J. M. Yeomans, Phys. Rev. 
E 59, 4366 (1999). 

[5] T. Ohta, H. Nozaki, and M. Doi, Phys. Lett. A 145, 304 
(1990); J. Chem. Phys. 93, 2664 (1991). 

[6] M. E. Cates, V. M. Kendon, P. Bladon, and J.-C. De- 
splat, Faraday Discuss. 112, 1 (1999). 

[7] N. P. Rapapa and A. J. Bray, Phys. Rev. Lett. 83, 3856 
(1999); F. Corberi, G. Gonnella, and A. Lamura, Phys. 
Rev. Lett. 81, 3852 (1998). 

[8] A. Coniglio and M. Zannetti, Europhys. Lett. 10, 575 
(1989); A. J. Bray and K. Humayun, Phys. Rev. Lett. 
68, 1559 (1992). 

[9] T. Ohta, D. Jasnow, and K. Kawasaki, Phys. Rev. Lett. 
49, 1223 (1982). 
[10] H. Orihara and Y. Ishibashi, J. Phys. Soc. Jpn. 55, 2151 
(1986). 

[11] S. M. Allen and J. W. Cahn, Acta Metall. 27, 1085 
(1979). 

[12] A. Cavagna and A. J. Bray, in preparation. 



4 



